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We study the phase diagram for a dilute Bardeen-Cooper-Schrieffer superfluid Fermi-Fermi mix- 
ture (of distinct mass) at zero temperature using energy densities for the superfluid fermions in one 
(ID), two (2D), and three (3D) dimensions. We also derive the dynamical time-dependent nonlinear 
Euler-Lagrange equation satisfied by the mixture in one dimension using this energy density. We 
obtain the linear stability conditions for the mixture in terms of fermion densities of the components 
and the interspecies Fermi-Fermi interaction. In equilibrium there are two possibilities. The first 
is that of a uniform mixture of the two components, the second is that of two pure phases of two 
components without any overlap between them. In addition, a mixed and a pure phase, impossible 
in ID and 2D, can be created in 3D. We also obtain the conditions under which the uniform mixture 
is stable from an energetic consideration. The same conditions are obtained from a modulational 
instability analysis of the dynamical equations in ID. Finally, the ID dynamical equations for the 
system are solved numerically and by variational approximation (VA) to study the bright solitons of 
the system for attractive interspecies Fermi- Fermi interaction in ID. The VA is found to yield good 
agreement to the numerical result for the density profile and chemical potential of the bright soli- 
tons. The bright solitons are demonstrated to be dynamically stable. The experimental realization 
of these Fermi-Fermi bright solitons seems possible with present setups. 

PACS numbers: 03.75.Ss, 64.75.-|-g, 03.75.Kk 



I. INTRODUCTION 



After the experimental realization of a trapped Bose- 
Einstein condensate (BEC) ,1,] there has been a great 
effort to trap and cool the Fermi atoms to degeneracy by 
sympathetic cooling in the presence of a second Bose or 
Fermi component. The second component is needed to 
facilitate evaporative cooling not possible due to lack of 
interaction in a single-component Fermi gas 0, Q . Apart 
from the observation of the degenerate Bose-Fermi mix- 
tures S'^Li 23Na-6Li and 87Rb-40K 0,I1, there 
have been studies of the following spin-polarized degen- 
erate Fermi-Fermi mixtures "'"K-'^^K and ^Li-^Li Q 
in different hyperfine states. 

Specially challenging has been the experimental real- 
ization of vortex lattice in a Bardeen-Cooper-Schrieffer 
(BCS) superfluid Fermi gas 0, M, HH, [11 in Bose-Fermi 
mixture employing a weak attractive interaction among 
the intra-species fermions by using a Feshbach resonance 
[isl . [l^ . This attractive interaction allows the formation 
of BCS pairs leading to a BCS superfluid [l^, [1^1 . In 
the last few years by further increasing this attraction 
several experimental groups have observed the crossover 
[l^ from the paired BCS state to the BEC of molecular 
dimers with ultra-cold two-hyperfine-component Fermi 
vapors of '^"K [13 and "^Li atoms [11 [H. Another possi- 
bility is to use two distinct Fermi atoms for this purpose 
as suggested in Ref . [l^l in a study of collapse in a Fermi- 
Fermi mixture (^Li-^^K is a possible candidate for future 
exploration.) The Feshbach-resonance management of 
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the Fermi interaction could be utilized to study a super- 
fluid Fermi- Fermi mixture in a controlled fashion [13[ . 

In Bose-Fermi mixtures, there have been several stud- 
ies on phase separation [2l|, [2^. 23l . [H, [2^ , soliton- like 
structures [IHl, and collapse [27 1. recently. The phase 
diagram of the Bose-Fermi mixture in three dimensions 
(3D) has been studied by Viverit et al. [1^, whereas 
the same in one dimension (ID) has been studied by 
Das 1241. Bright solitons have been observed in BECs 
of Li |28l| and Rb [2§| atoms and studied subsequently 
[30| . It has been demonstrated using microscopic [sj 
and mean-field hydrodynamic [1^ models that the for- 
mation of stable fermionic bright solitons is possible in 
a degenerate Bose-Fermi mixture in the presence of a 
sufficiently attractive interspecies interaction which can 
overcome the Pauli-blocking repulsion among fermions. 
The formation of a soliton in these cases is related to the 
fact that the system can lower its energy by forming high 
density region (bright soliton) when the interspecies at- 
traction is large enough to overcome the Pauli-blocking 
interaction in the degenerate Fermi gas (and any possible 
repulsion in the BEC) [32|. There have also been studies 
of mixing-dcmixing transition in degenerate Bose-Fermi 
[sst and Fermi-Fermi (sj mixtures, and soliton formation 
in Fermi- Fermi mixtures [35| . 

In this paper we investigate the phase diagram of a 
BCS superfluid Fermi-Fermi mixture of fermion compo- 
nents of distinct mass at zero temperature using energy 
densities for the superfluid Fermi components in one, two 
(2D), and three dimensions. We derive the conditions of 
stability of the mixture in terms of the densities of the 
components and the strength of interspecies interaction. 
The two possible phases of the mixture are a uniformly 
mixed configuration and a totally separated pu re-p hase 
configuration. Unlike in a Bose-Fermi mixture [231 . [2^ , 
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no complicated mixed phases are allowed in a superfluid 
Fermi-Fermi mixture in ID and 2D. However, a mixed 
and a pure phase is allowed in 3D. In ID, two pure and 
separated phases of the fermion components appear for 
low fermion densities, whereas the opposite is found in 
3D. In ID, a uniform mixture appears for large fermion 
densities with the opposite taking place in 3D. In 2D, 
the condition for uniform mixture and phase separation 
is independent of density of the components. In ID, we 
find the uniform mixture to be unstable for small fermion 
densities, whereas In 3D, the uniform mixture is unstable 
for large fermion densities. 

The ID configuration is of special interest due to soli- 
ton formation by modulational instability of a uniform 
mixture. To study this phenomenon we derive a set of 
dynamical equations of the system as the Euler-Lagrange 
equation of an appropriate Lagrangian. The condition 
of stability of the uniform mixture and the formation of 
soliton for attractive interspecies Fermi- Fermi interaction 
were studied from an energetic consideration as well as 
with a linear stability analysis of the constant-amplitude 
solution of the above dynamical equations. We solved 
the ID dynamical equations numerically and variation- 
ally to study some features of the bright solitons. The 
numerical results for the density of the fermion compo- 
nents as well as their chemical potentials are found to be 
in good agreement with the variational findings. These 
bright solitons are found to be stable numerically when 
they are subjected to a perturbation. 

The dependence of Fermi energy densities in ID and 
2D on atomic densities has counterparts in Bose sys- 
tems and the analysis presented here is also applicable 
to these Bose systems. The 2D Fermi energy density has 
a quadratic dependence on atomic density as in a dilute 
BEC obeying the Gross-Pitaevskii equation, thus allow- 
ing the present results to be applicable to such a BEC 
The ID Fermi energy density, on the other hand, 
has a cubic dependence on atomic density as in a Tonks- 
Girardeau jjg] (TG) Bose gas observed recently [37[, thus 
making the present results applicable to this system. 

The paper is organized as follows. In Sec. II we con- 
sider the stability condition of a uniform BCS superfluid 
Fermi-Fermi mixture from an energetic consideration. In 
Sec. Ill we consider a two-phase BCS superfluid Fermi- 
Fermi mixture in ID, 2D and 3D and study the possi- 
bility of the formation of two phases from a consider- 
ation of pressure, energy and chemical potential of the 
system. We can have two pure phases or a uniformly 
mixed phase in all dimensions. In addition, in 3D, we 
can have a pure and a mixed phase. In Sec. IV we con- 
sider the Euler-Lagrange nonlinear dynamical equations 
for the system in ID and study the modulational insta- 
bility of the constant-amplitude solution representing the 
uniform mixture. The condition of modulational insta- 
bility for attractive Fermi-Fermi interaction is found to 
be consistent with the condition of stability of the uni- 
form mixture obtained from an energetic consideration in 
Sec. II. We further solve these dynamical equations nu- 



merically and variationally to analyze the properties of 
the Fermi-Fermi solitons. Finally, in Sec. V we present 
a summary of our study. 



II. UNIFORM SUPERFLUID FERMI-FERMI 
MIXTURE 

A. Energy Density of a Component 

We consider a single-component dilute BCS superfluid 
of spin-half Fermi atoms of mass m and density with 
a weak attraction between fermions with opposite spin 
orientations. In 3D, the energy density of this system is 
given by [H, [H, H, HI , 
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where ep = (fikp) /(2m) is the Fermi energy, fikp 
is the Fermi momentum (this expression was first ob- 
tained by Lee and Yang [33 in the weak-coupling BCS 
limit). Modifications to this expression for a description 
of the BCS-BEC crossover, for stronger attraction be- 
tween fermions, have also been considered [4l] |. The total 
density of the fermions in a 3D box is obtained by filling 
the quantum states up to the Fermi energy and is given 

by 713 = 2(2^)-3/o'=^47rfc2dfc = (3^2)-i {2meF/n^f'''. 
(The factor of 2 in the expression for ns accounts for 
BCS pairing in each level.) Hence the energy density in 
([T|) becomes 
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with .43 = ?i2(37r2)2/3/(2m). 

Similarly, the energy density of a dilute ID superfluid 
of atom density ni is given by (4^ . l43j 
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This was obtained using the Gaudin-Yang (GY) model 
[4^ of fermions weakly interacting via zero-range {5- 
function) potential, and was later extended to the de- 
scription of the BCS-to-unitarity crossover [i^]. (For re- 
pulsive interaction the GY model gives [i^l a Tomonaga- 
Luttinger liquid [Ssj . while for attractive interaction it 
leads to a Luther-Emery liquid [i^. For weak attrac- 
tion the ground state of the system is a BCS superfluid 
p^ . [47} . With the increase of attraction, the strong- 
coupling regime of tightly bound dimers is attained, 
which behaves like a hard core Bose gas, or like a ID non- 
interacting Fermi gas, known as the TG gas [3^,|4^ .) The 
general solution for the ground-state energy in the GY 
model has been obtained by solving the Bethe ansatz [i^ 
equations for all strengths of 5 interaction connecting the 
weak-attraction regime of BCS condensate to the strong- 
attraction regime of of tightly bound dimers described 
by the Lieb-Liniger model [50j of repulsive bosons. This 
solution can be presented as an expansion series in limits 
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of weak or strong interactions. The limiting value of this 
solution in the weak interaction BCS limit is given by 
Eq. © pin. 

The fermion density of the BCS superfluid in a ID box 
is ni = 2(27r)~^ I-kp = (2/7r?i) y^2meF, hence, in this 
case, ep = n^h'^nl/{8m), and energy density ^ becomes 
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with Ai = h^TTy{8m). The energ y d ensity of a TG gas 
(36| is given by Stg = ?i^7r^nf /6m j48| and is very similar 
to that given by Eq. The difference in numerical 

factors between the two expressions is due to pairing in 
the present Fermi superfluid allowing two fermions (spin 
up and down) in the same quantum level. Hence the ID 
results of the present study is also applicable to a TG 
gas. 

Finally, a counterpart of relations ([T]) and (|4]) for the 
2D superfluid is [52| £2D = (l/2)n2eF, the 2D den- 
sity being n2 = 2(27r)~^ jj"^ 2Trkdk = (m/iih?^ ep, with 
Sp = Trh^n2/m. Thus, the energy density of the 2D su- 
perfluid can be written as (5^ 
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T'2 = 2^2n2, 
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with A2 ~ nh^ /m. 

Here we specify the criteria of applicability of Eqs. ([2]), 
Q, and (O for different dimensionalities. These results 
are valid for a dilute BCS superfluid. In 3D, at low den- 
sities, fci?|ai?| << 1 with ap the Fermi-Fermi scattering 
length, gaps are small and have little effect on the total 
energy of the system [4^ . The total energy density of the 
ground state can then be expanded in powers of the small 
parameter kp\ap\. At low densities Eg. ^ includes the 
lowest order term in this expansion [39[. The condition 
fcF|a_F| << 1 of validity of Eq. ^ can be related to 
the gas parameter nsla^p in 3D: n3|a_Fp << l/(37r^), 
as the density 713 = fc|,/(37r^). In ID, for a 6 interac- 
tion of strength gi the dimensionless coupling constant 
7 = mgi / {h^ ni) and the condition of validity of Eq. ^ 
is I7I << 1 [4^ . In two dimensions an attractive interac- 
tion leads to a bound state of energy eo and the condition 
of dilutencss for the validity of Eq. ^ can be expressed 
as eo/ep « 1. 



B. Stability Condition of the Uniform Mixture 

We consider a uniform mixture of two types of 
fermions, containing iV^, z = 1,2, atoms (of mass mi = m 
and 1712 = m/X) , in a box of size S (in ID the size is 
a length, in 2D an area, and in 3D a volume) with dis- 
tinct mass at zero temperature. The energy density of 
the uniform mixture is given by 
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respectively, for ID, 2D and 3D systems, where nj^^i-j = 
Ni/ S denotes the density of each component in dD, d — 
1,2,3. The nonlinear terms involving 1712 = 47r 71^012/77112 
in above equations represent the interaction between two 
types of atoms arising solely from the atomic scattering 
length ai2, where 77112 is the reduced mass of atoms. The 
terms involving Aj in the above equations, although is 
similar to the gn^ jl interaction term for bosons (with 
g ~ A-nh'^a/m representing the self interaction of a di- 
lute boson gas with a the Bose-Bose scattering length 
and m the mass of an atom), have a different origin as 
we have seen. These terms originating from the energy of 
the fermions occupying the lowest quantum levels at zero 
temperature obeying Pauli principle generate an effec- 
tive repulsion between the fermions and is usually called 
Pauli-blocking interaction. 

The chemical potentials = d£/dni for species i ~ 
1,2 in ID, 2D and 3D, are given, respectively, by 

fii = ^l?^?(l) +5i2"-i(2), M2 = .9i2n-i(i) + AAin^2), 

(9) 

Ml = -42ri2(l) + 5l2'T'2(2), A'2 = .9l2'T-2(l) + '^^2"2(2), 

(10) 



Ail 



= A 
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3"3(r) + 5l2"-3(2), A*2 = 5l2'T-3(l) + -^^3^3(2) • 
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The uniformly mixed phase is energetically stable if 
its energy is a minimum with respect to small variations 
of the densities, while the total number of fermions and 
bosons are held fixed. The conditions of stability [are the 
conditions of a minimum of £(77(1), 77(2)) as a function of 
two variables n(^i-^ and ?7(2) and] are given by 
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>0, 
(13) 



where we have dropped the dimension suffix. The solu- 
tion of these inequalities gives the region in the parame- 
ters' space where the uniformly mixed phase is energet- 
ically stable. Using Eqs. (|12|) and (fT3ll the condition of 
stability of the uniform mixture in ID, 2D, and 3D are 
given, respectively, by [23l.[4l| 



4AiA77(i)77(2 

AlX>gj2, 



AAiX > 9gi2n' 
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(14) 
(15) 
(16) 



These conditions are determined by gi2 and not the sign 
of .gi2. 
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In ID, we find from Eq. lfT4|) with a finite j that at 
small fermionic densities (small n(i) and n(2)) the uni- 
form mixture is unstable: the ground-state of the system 
displays demixing if 912 > and becomes a localized 
Fermi- Fermi bright soliton if 1712 < [i^. The mixture 
is stable at large fermionic densities. In 2D, Eq. (fT5|) 
reveals that the condition for stability is independent of 
density. In 3D, Eq. (|16p predicts that for a finite g^2, the 
mixture is unstable at large fermionic densities, leading 
to collapse for gi2 < and to demixing for gi2 > 0, and 
stable at small fermionic densities. It is realized that as 
we move from ID to 3D through 2D, the condition of sta- 
bility of the uniform mixture changes from large fermion 
densities to small fermion densities. This result is quite 
similar to that in a Bose- Fermi mixture [l^, , where 
the condition of stability of the uniform mixture is inde- 
pendent of the bosonic density and has a similar depen- 
dence on fermion density, e.g., during the passage from 
ID to 3D through 2D, the condition of stability changes 
from large fermion density to small fermion density. 

From inequality the stability condition of a single 
component uniform gas can be represented as dfii / dni > 
0, which, using Eqs. pU)) and PT|) . is reahzed for 

Ad > denoting a repulsive system. In the presence of a 
second component, inequality can be written as [2^ 

|^-f|^)'^>0, (17) 

as d^2/drL(i) = 9^i/9n(2). The first term d^i/drni^ in 
inequality pT]) represents the effective repulsion among 
fermions of type 1. The second term, representing an 
induced interaction due to the presence of component 2, 
reduces the repulsion and tries to destabilize the uniform 
mixture. The uniform mixture becomes unstable when 
the second term in inequality (jl7p becomes larger than 
the first term. This happens for both attractive and re- 
pulsive interspecies interaction gi2 = 9/i2/9n(i). 

The inequality (jl3p can be written as 

clcl> 4:gj2n (^1)71(^2), (18) 

where Ci = ^2ri(j)(9^i/9n(j)) represent sound velocities 
in the two superfluid components, i = 1,2. The sound 
velocity C12 of the ID Fermi-Fermi mixture can be ob- 
tained following a procedure suggested by Alexandrov 
and Kabanov [4^, |53l for a two-component BEC: 



Cl2 = -^]J 4 + 4 ± ^ (cf - c2)2 + 165j2'^(l)"(2) • 

(19) 

The homogeneous mixture becomes unstable when the 
sound velocity C12 becomes imaginary, e.g., when inequal- 
ity is violated. 



III. TWO-PHASE SUPERFLUID FERMI-FERMI 
MIXTURE 

In the last section we considered a uniform mixture 
of two components in equilibrium. Here we explore the 
more interesting case of two types of fermions with dif- 
ferent possible densities in different regions of a box of 
size S. The components may mix uniformly or form sep- 
arate phases depending on the initial conditions — mass, 
density, interspecies interaction etc. 

The conservation of the number of particles, Ni and 
A^2, of the two species can be expressed as [2^ IIJ] 

2 2 
= Sn, = 5^ = 1' (20) 

where i = 1 , 2 represent the species and j = 1,2 repre- 
sent the phases (different region with distinct density of 
gas), rii = Ni/S represent the overall density of the two 
species, riij is the density of species i in phase j, and 
Sj = Sfj represent the size of each phase with fj the 
fraction of size in phase j. For a two-component system 
one can have only two distinct phases, j = 1, 2, as the 
inclusion of more phases leads to inconsistency [23| . Here 
we have dropped the dimension label d and also removed 
the parentheses () from the component label i. 
The total energy of the system is given by 

2 2 

£; = ^£;,^5^/,£:, (21) 

where Sj denotes the energy density of phase j and Ej 
its total energy. The pressure Pj of phase j is given by 
Pj ~ —dEj/dSj. The chemical potential of component i 
in phase j is defined by fj,ij = d£j / driij. 

For equilibrium, the pressure in one phase has to be 
equal to that in the other. If two phases are occupied by 
atoms of the same type, the chemical potential for that 
type of atoms in two phases should also be equal so that 
the equilibrium can be energetically maintained. If the 
atom density of one type of atom in a phase is zero then 
the chemical potential of that type of atom in this phase 
should be larger than that in the other, so that the atoms 
do not flow to the phase with no atoms of this type [23j . 

In the following we consider a system composed of two 
phases comprising of fractions /i = / and /2 = (1 — /) 
of size S. There are three following possibilities to be 
analyzed in ID, 2D, and 3D, although some of them may 
not materialize in a particular case: 

(i) Two pure and separated phases with one type of 
atom occupying a distinct phase. 

(ii) A mixed and a pure phase where the density of one 
type of atoms is zero in one phase. 

(iii) Two mixed phases where both phases are occupied 
by both type of atoms. 

In the following we deal with the three possibilities in 
ID, 2D and 3D. First, we consider the 2D case as the 
algebra is significantly simpler in this case. 
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A. Two-Dimensional (2D) Mixture 

From Eqs. ([7]) we find that the expressions for total 
energy and pressure in this case are 



SjSj — Sj 



In deriving Eq. we recall that 



l/S-j. From 



Eq. (p2|) the chemical potentials are given by 

/il J = ^2711 J + 312^2^ (24) 



M2j = ^2^712 J + 5l2"l J 



(25) 



where we have used Eqs. ([30]) . The energy of the phase- 
separated system with the same number of atoms is 



E, 



scp 



s 



+ ]-A2\{\ - f)nl 



(33) 



Using Eq. (|26p . one has for the difference 

Smix - Sscp = Sf{\ - /)n2 2yA[5i2 - ^2%/A]. (34) 

When i?mix > ^-scp the system naturally moves to the 
separated phase and this happens for > ^i-^i consis- 
tent with inequality leading to a stable separated 
phase. In the opposite limit, when i?mix < ^-sop, the en- 
ergetic consideration favors the uniform mixture and this 
happens for < A\\, consistent with inequality ([15]). 



1. Two Pure Phases 

In case of two pure and separated phases one should 
have, for example rii 2 — 712.1 = corresponding to the 
type one atoms occupying phase 1 only (rii^i 7^ 0) and 
type 2 atoms occupying phase 2 only (712,2 7^ 0). 

Equality of pressure Pi = P2 in the two phases yields 



'1,1 



A772 2- 



(26) 



As the number of atoms is zero in one of the phases, 
one has the inequalities ^2,2 < M2.1 a-nd /ii,i < ^1.2 on 
the chemical potential, which, using Eqs. (|24p and (P5)l . 
become 



^2An2,2 < 5l2"l,l, 
^2711, 1 < 512712, 2- 



(27) 
(28) 



Eliminating ni.i and 712.2 among Eqs. (|26p . (|27p . and 
([28)1 we get 



(29) 



consistent with inequality (|15p . We have the uniform 
mixture for inequality (|15p : for the opposite inequality 
P9|l we have the separated phases in equilibrium. These 
inequalities are independent of the atomic densities. 

In the present case the overall densities of the two 
species are given by 

"1 = /l"-l,l = /"-1,1, «2 = /2n2,2 = (1 - /)"-2,2- (30) 

Let us now consolidate these findings using energetic 
considerations comparing the total energy of a phase- 
separated configuration with that of a uniform mixture. 
The energy of the mixture is given by 
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2. A mixed and a pure phase 

Here we consider one mixed phase (phase 1) and one 
pure phase (phase 2) consistent with 711^2 = 0, which 
means that the type 1 atoms occupy only phase 1, 
whereas type 2 atoms occupy both phases 1 and 2. Using 
Eq. ([23]) the equality of pressure in two phases leads to 



]^A2n\-^ + 5l2lll,lll2,l + ^^2A7l2^i = ^^2A7l2^2- (35) 

From Eq. (|25p the equality of the chemical potential of 
type 2 atoms in two phases (/i2,2 = ^2,1) leads to 



711,1 = ^2A(?l2,2 - ll2,l)/ff 



12- 



(36) 



From Eq. (|24|) the inequality of the chemical potential of 
type 1 atoms in two phases (/ii.i < /^i,2) leads to 

^2711,1 < 512(12,2 - 712, 1), (37) 

which using Eq. (|36p yields 

AlX < gl2. (38) 

Substituting Eq. ((36)) into Eq. ((35l) and after some 
straightforward algebra we obtain 



(A2A-522)(n2,i- 712,2)^ = 0, 



(39) 



which allows two possibilities. For A\X ^ 5J2, the only 
solution is the trivial, nevertheless unacceptable, one 
7i2,i = 7i2,2, which means that the type 2 atoms form 
a uniform configuration and not a mixed phase. How- 
ever, if AjA = 5^2! O'^'^ <3an. have a mixed phase with 
712.1 7^ 7i2,2- Nevertheless, this condition enters in con- 
tradiction with inequality psp , showing that one cannot 
have one mixed and one pure phase in this case. 

Next we consider the possibility of two mixed phases. 
The equality of pressure and chemical potential of each 
species in two phases leads to the following conditions 



1 



1 



7;Mnii~<2) + o^2A(7i^,i 



= 512(111, 2"'2, 2 - 111, lll2,l), (40) 
^27ll,l -1-51212,1 = ^2711,2+5121-2,2, (41) 
^2A7l2,l -|-5l2ll,l = ^2A7l2,2 + 5l2ll,2- (42) 



This set of equations have only the trivial solutions ni^i = 
ni_2 and 712.1 = ^^2,2 corresponding to uniform mixture. 
Hence two mixed phases cannot be in equilibrium. 



B. One-Dimensional (ID) Mixture 

From Eqs. ([6]), wc find that the expression for total 
energy and pressure in this case are 



dEj 2 n 2 o 



From Eq. the chemical potentials are given by 



Mij = + .9l2"-2,i 

M2j = AiXnl + 312^1, 



(43) 

r 
(44) 

(45) 
(46) 



1. Two pure phases 

In case of two pure and separated phases one should 
have, for example, n2,i = for phase 1 and ni 2 = for 
phase 2. The condition of equal pressure then yields 



Ann 



(47) 



For equal-mass fermions A = 1, and one obviously have 
have the trivial solution ni 1 = 712.2 or the densities of the 
two species are equal. Of course, for A 7^ 1 the densities 
of the two species could be different. Chemical potential 
condition /i2,2 l£ M2.1 yields 

AiXnl ^ < .9i,2rii,i- (48) 

Chemical potential condition fii i < 111^2 yields 

Ain\ 1 < .gi,2n2,2- (49) 



Eliminating Tii.i between Eqs. (|47p and or between 
Eqs. (gT]) and ^ we get 



n2,2<Bi, Si =.gi2/(^iA2/3). 



(50) 



From Eqs. (|T7)) and ([50)1 we obtain the following restric- 
tion on 71 1 1 



711 



,i<Cu Ci =,9i2/(AiAi/3). 



(51) 



In this case a phase diagram showing the total densi- 
ties of type 1 and 2 fermions for which the system can 
completely separate, can be obtained from Eq. ([5(7|l if 
we allow / to vary from to 1 and use conditions ([501) 
and ([5T|) . This is illustrated in Fig. [TJ The light gray 
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FIG. 1: (Color online) Phase diagram for Fermi-Fermi mix- 
ture in one dimension (ID). The plotted density ni is in units 
of Ci = 5i2/(AiA^/^) and n2 in units of Bi = giiKAiX^''"^). 
The light gray area represents two pure and separated phases 
while the dark gray area represents the domain of stable uni- 
form mixture. Below the curved line in the clear area the 
uniform mixture is unstable against small fluctuations. 



area represents pure phases and the dark gray area repre- 
sents the stable uniform mixture. The uniform mixture 
is unstable in the clear area below the curve given by 
inequality (HH). For attractive interaction, one has the 
formation of bright solitons by modulational instability, 
(discussed in Sec. IV). For repulsive interaction one can 
have a partially demixcd configuration in the clear region 
in Fig. □ 

Now let us see if the system spontaneously move into 
the phase-separated configuration from an energetic con- 
sideration. The energy of the mixed system is 



E„ 



S 



S 



1 o 1-1 

-AlTlj -I- 711712312 + 2 AAl7l2 



Ainl + 71i,i7l2,2./(l - /)gi,2 



-A,nl,il-f)'X 



(52) 



(53) 



Equation ([SS]) is obtained with the use of Eq. pO|) . The 
energy of the separated phase system with the same num- 
ber of atoms is 

Esep = lA,S[nlJ + nl^Xil - /)]• (54) 

Using Eq. (|47p . and after some straightforward algebra, 
the difference A = (-©mix ~ Egcp) is given by 

A = 7l2 5Al/3[gi2 - Ai7l2,2A2/3]/(l - /) (55) 
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Considering the restrietion ((50)) in the separated phase, 
Eq. ([55]) yields the foUowing inequahty 

A = nlJil - /)5Ai/3[i _ n2,2/Bi]gi2 > 0. (56) 

For density ranges where equihbrium is possible / 7^ 
and / 7^ 1, -Bsop is always less than i?mix- Hence, ener- 
getically the two species of fermions can separate. 



This set of equations have only the trivial solutions ni^i = 
2 and 712^1 = n2_2 corresponding to uniform mixture 
and that is also possible when the condition of uniform 
mixture is satisfied. Hence two mixed phases cannot be 
in equilibrium. 

C. Three-dimensional (3D) Mixture 



2. A mixed and a pure phase 

Now let us consider a mixed phase (phase 1) and a 
pure phase (phase 2) and consider the case ni^2 ~ 0. 
The equality of pressure now leads to 

2 2 2 

-Ainl ^ + -AiXnl -^ + gi2ni,in2,i = -^lAn^ 2- (57) 

The equality of chemical potential of species 2 in two 
phases (/i2,i = ^2,2) yields 



711,1 = AiX{n\ 2 - n\ i)/gi2- 



(58) 



Eliminating ni.i between equations (|57p and (|58p (after 
some straightforward algebra) we get 



(59) 



where x = n2.i/n2.2, A = {712,2/ . After cancelling 
the trivial factor (1 — x^ from both sides of Eq. ([59|l . we 
get 



2A(l + a;)3(l-x) = x + 2 



(60) 

From Eq. ([60)) we find that the solution a; = is obtained 
for A = 1 corresponding to 7^2,1 = ni,2 = 0, 712,2 = Bi 
and 711,1 = C'l. The densities of the first component are 
Til, 2 = and 7ii,i = Ci. This is the special case con- 
sidered in Sec. IHBl [see, Eqs. §0^ and §1^]. The 
solution 7i2,2 = -Bi (A = 1) is a solution of two pure 
phases corresponding to a; = 0. The domain of solution 
of mixed phase corresponds to 712,2 > -Bi (A > 1) cor- 
responding to X > (recall that the fraction x cannot 
be negative.) Hence for the present mixed phase to exist 
Eq. ([50]) should have the solution x +0 for A ^ +1. 
However, we find from Eq. ([60]) as A is made slightly 
greater than 1, the solution a; = turns negative (un- 
physical). [Please note that for A = 1 Eq. (|5D|) has two 
real roots: a; = 0, 0.7399...; the latter (spurious) root 
is not of present physical interest.] Hence we conclude 
that a mixed and a pure phase cannot be realized in the 
present mixture. 

Finally, one can consider the possibility of two mixed 
phases. The equality of pressure and chemical potential 
of each species in two phases leads to 



3^i«i -"L) 



^l'^l,l + 5l2'^2,l 

AXnl 1 + 312111, 1 



-AiXinl, 



2) 



.9l2("-l,2'T'2,2 - "l,l"2,l), (61) 
,2 -T-5l2"-2,2, (62) 



Ainf. 



AXn2^2 + 5i2"-i,2- 



(63) 



From Eqs. ([8]), we find that the expression for total 
energy and pressure in this case are 



Bj — SjSj — Sj 



/I 5/3 I I ^ ,1 \ 5/3 



Pi 



dEj 2 5/3^ ^2 5/3 



From Eq. ([22]) the chemical potentials are given by 



^l,j = ^3%f/ + 5l2«2,i 

2 /3 

M2,j = ^3An2j- +gi2iii 



(64) 
(65) 

(66) 
(67) 



1. Two pure phases 

Again for two pure and separated phases we take 
Til 2 = 7i2,i = 0. The condition of equal pressure in two 
phases then leads to 



5/3 _ ^ 5/3 

^1,1 ~ ^'^2,2 ■ 



(68) 



The chemical potential condition ^2,2 < M2,i yields 

"-1,1 > AA3n2^2^/gi2. (69) 
The chemical potential condition /ii i < ^1,2 yields 

"2,2 > A3?lJ/f /ffl2. (70) 

Eliminating iii 1 between Eqs. ([55]) and (p5)) or between 
Eqs. (inSl) and dTD]) we obtain 



,2 > S3, B3 = (A2/5A3/ffi2)' 



(71) 



«2 

Similarly, eliminating 112,2 between Eqs. (j68p and ((691 
we get 



11,1 > C3, C3 = (A3/5^3/gi2) = 



(72) 



In this case a phase diagram showing the total densi- 
ties of type 1 and 2 fermions for which the system can 
completely separate, can be obtained from Eq. ||5D|) if we 
allow / to vary from to 1 and use conditions (|7T|) and 
dZH). This is illustrated in Fig. [H 

To see the separation of the two types of fermions from 
an energetic consideration, we calculate the energies of 
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For 1 > / > 0, the quantity given by (|77|) is always 
positive. Hence the separated phase has less energy than 
the mixed phase and the system will spontaneously move 
into the phase separated configuration. 

In this case also two mixed phases cannot be in equi- 
librium as in ID. 



2. A mixed and a pure phase 

Again we consider a mixed (species 2) and a pure 
(species 1) phase and consider the case ni^2 = 0. The 
equality of pressure now leads to 

^ /I 5/3 I I ^ /I \ 5/3 2 5/3 

-Asn^'j + 312^11, 1*^2,1 + g^3Ari2,i = g^3'^"-2,2 • (78) 

The equality of chemical potential of species 2 in two 
phases (/i2,i = ^^2,2) yields 



FIG. 2: (Color online) Phase diagram for Fermi-Fermi mix- 
ture in three dimensions (3D). The plotted density ni is 
in units of C3 = (A^^^As/gig)^ and n2 in units of Bs = 
(A^/^As/gig)'^. The light gray area represents two pure and 
separated phases while the dark gray area represents the do- 
main of stable uniform mixture. The clear area represents 
the region where a mixed and a pure phase can exist. Above 
the curved line the uniform mixture is unstable against small 
fluctuations. 



the mixed and separated configurations. The energy of 
the mixed phase is [2^ 



E„ 



S 



= S 



5/3 
2 



(73) 



-Ag/^/^n^/f +512^1,1^2,2/(1-/) 



5-. -2,2' 

The energy of the separated phase is 



3 . 5/3 J. , 3 . 5/3., ^^ 



(74) 



(75) 



Using Eq. ((68|) the difference A = (fmix — ^scp) can be 
written as 

A = 5[3A3An^f|/(/'/' - l)/5 + .9i2nLA3/5/(l - /) 
+ 3AsX4^^{l - /){(! - ff/' - l}/5]. (76) 

Using inequality (|70|) . Eq. (|76|) yields 



A > ^ 



9l2 



-J{1 f") 
5 



/(I -/)- 5(1 -/){1-(1 -/)'/'} 



(77) 



ni,i = yl3A(n: 



2/3 
2,2 



2/3 



t)l9i 



(79) 



Eliminating rii i between Eqs. ([75)) and ([75]) and after 
some straightforward algebra we get 

2A(1 + xf/^ - (1 - xY'^{Zx^ + Qx^ + 4.T + 2) = 0,(80) 

where x — (^2, 1/712,2)^^'^, and A = (n2,2/i?3)~^/^. From 
Eq. ([5D|) we find that the solution a; = is obtained 
for A = 1 corresponding to 71,2,1 = 0, 772,2 = -S3, '^i,i = 
^3, 77i,2 = 0. This is the limiting case of two pure and 
separated phase studied in Sec. IIICl. (In addition for 
A = 1, Eq. has the spurious or unphysical root 

x — 0.90278..., which we do not consider here.) For two 
purely separated phases we have seen that 772,2 > S3 
whence A = {112,2/ B^)^^^^ < 1- The domain for a mixed 
and a separated phase then should have A > 1. To find 
this domain we solve Eq. ([50]) for x > using different 
A. Such solutions appear in the range 1.217... > A > 1. 
Using this solution for x we obtain 772,2 and 772,1 from 
the definitions of A and x, respectively. Finally, 771,1 is 
obtained from Eq. ((7^ . The results so-obtained for 772,2, 
772,1, and 77 1,1 for different A are used in 



771 = /77l,l, and 7i2 = /772,1 + (1 - /)772,; 



(81) 



to calculate the domain of 771 and 772, by varying / in the 
range 1 > / > 0, which allows a pure and a mixed phase. 

We show the 3D phase diagram for total densities of 
type 1 and 2 fermions in Fig. (5] In this figure the light 
gray area represents the domain of two separated phases 
and the clear area that of a mixed and a separated phase 
as calculated above. The remaining dark gray area repre- 
sents the domain of stable uniform mixture. The uniform 
mixture is unstable above the curve given by Eq. (|16|) . 
Qualitatively, Fig. [5] is quite similar to Fig. 3 of Vivcrit 
et al. [2^ for a Bosc-Fcrmi mixture. 

If wc compare Figs. [T]and[2]wc find that in ID the pure 
phases appear at small densities, and uniform mixture at 
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large densities. The uniform mixture is stable at larger 
densities. The opposite happens in 3D. If we compare 
the findings of Viverit et al. [23| for a study of the phase 
diagram of a Bose- Fermi mixture in 3D and compare with 
the study of Das in ID we find that such an inversion 
also takes place there. Moreover in ID there cannot be 
a mixed and a pure phase for a Fermi-Fermi mixture, 
which is possible in 3D. 



IV. DYNAMICAL EQUATIONS IN QUASI-ID 
SUPERFLUID FERMI-FERMI MIXTURE 

A. The Model 

Of the three dimensional possibilities — ID, 2D and 3D 
— the ID case deserves special attention. In ID, if the 
interspecies Fermi-Fermi interaction is attractive, in the 
domain of instability of the uniform mixture one can have 
the formation of bright soliton by modulational instabil- 
ity. To perform a careful study of the nature of these 
bright solitons (and their dynamical stability) we derive 
the Euler-Lagrange equations in ID from its Lagrangian 
density. 

We consider a mixture of A^i superfiuid atomic 
fermions of mass toi(= m) and N2 superfiuid atomic 
fermions of mass m2(= m/A) at zero temperature 
trapped by a tight cylindrically symmetric harmonic po- 
tential of frequency ujj_ in the transverse (radial cylin- 
dric) direction. We assume factorization of the trans- 
verse degrees of freedom. This is justified in ID con- 
finement where, regardless of the longitudinal behavior 
or statistics, the transverse spatial profile is that of the 
single-particle ground-state [3, [s^, HH] • The transverse 
width of the atom distribution is given by the character- 
istic harmonic length of the single-particle ground-state: 
a±j = ^h/ (rnjijj±), with j = 1,2. The atoms have an ef- 
fective ID behavior at zero temperature if their chemical 
potentials are much smaller than the transverse energy 
hLLi± [23, [13, \EM ■ The interspecies Fermi-Fermi interac- 
tion is characterized by a contact potential with scatter- 
ing length ai2, which can be repulsive or attractive. 

We use a mean-field Lagrangian to study the static and 
collective properties of the ID superfiuid Fermi-Fermi 
mixture as in the Ginzburg-Landau theory [T^ . The La- 
grangian density £ of the mixture reads 



(82) 



The term is the fermionic Lagrangian for component 
i, defined as 



' dt 



2m 



dipt 



dz 



^1* 



local probability density of the ith component. Here m 



is the effective mass of superfiuid flow in the Ginzburg- 
Landau theory. There is experimental evidence that 
this effective mass is twice the fermion mass (to^^ = 2mi) 
and we shall use this effective mass in the following study. 

Finally, the Lagrangian density £12 of the interaction 
between the two Fermi components is taken to be of the 
following standard zero-range form (ssl . [53 | 

£12 = -5i2|^i|V2|' , (84) 

where 512 = 2TiLo^ai2 is the ID Fermi-Fcrmi interaction 
strength. 

The Euler-Lagrange equations of the Lagrangian £ are 
the two following coupled partial differential equations: 



ihdti'i 
iWti'2 



Ami 



4m2 



.912 ?^2 



5.^ + 4 «2+. 912^^1 



i'2 



(85) 
(86) 



with the normalization ji/'ipdz = Ni. 

It is convenient to work in terms of dimensionless 
variables defined in terms of a frequency lj and length 
I = y/Ti/ {2miuj) by ipj = V^j/v^: ^ = 2t/tj, z = zl, and 
512 = 5i2?i^/(4mi^). With these new variables Eqs. ([85]) 
and (18611 can be written as 



idtipi = 



idtil^2 



-9 



.9l2"2 



Xdl + XAnl + .91271.1 



(87) 



V^2 , (88) 



where A = and where we have dropped the hats 

over the variables, and where A = 7711/7712, ni = \ipi\'^, i = 
1,2 with the normalization J^^n.idz ~ Ni. Equations 
([57)1 and with diagonal quintic nonlinearity are the 
equations satisfied by two coupled TG Bose gas [s^ and 
hence the analysis of Sec. IV also applies to a TG gas. 

For stationary states the solution of Eqs. ([57)1 and fSB]) 
have the form ipi = (f>i exjp{—ifiit) where fXi are the respec- 
tive chemical potentials. Consequently, these equations 
reduce to 



Aii^i 



M202 



- 9^ + Anl + 5i2?i2 



Aa2 



XAnl 



5l2"-l 



(89) 



(90) 



A repulsive interspecies Fermi-Fermi interaction is pro- 
duced by a positive gi2, while an attractive Fermi-Fcrmi 
interaction corresponds to a negative gi2. 



(83) 



where A^*"* = h'^n'^ / {8mi), tpi^z^t) is the field of the ith 
component of the BCS Fermi superfiuid along the lon- 
gitudinal axis, such that ni(z,t) = \'ipi(z,t)\'^ is the ID 



B. Modulational Instability 

To study analytically the modulational instability (35l . 
[56| of Eqs. ([57)1 and ([55)) we consider the special case of 
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attractive Fermi-Fcrmi interaction while these equations 
reduce to 

idti^i = +A|V'i|"-5i2|V'2|']^i , (91) 

^^t^P2 = [-Xdf^+AX\iP2f-giM'']^2, (92) 

where we have taken the interspecies interaction to be 
attractive by inserting an expHcit negative sign in gi2. 

We analyze the niodulational instability of a constant- 
amplitude solution corresponding to a uniform mixture in 
coupled Eqs. (PT|) and ([5^ by considering the solutions 

ipio = Hioexp(i5i) EE Hioe**(^'i^^^o-'^^io\ (93) 



<P20 



H2oexp(zJ2) = H2oe**(^^^'''"-^'^''™^ 



(94) 



of Eqs. (PTj) and ([5^ . respectively, where H^o is the am- 
plitude and Si a phase for component i. The constant- 
amplitude solutions, describing an uniform mixture, de- 
velop an amplitude-dependent phase on time evolution. 
We consider a small perturbation HiCxp(i(Si) to these so- 
lutions via 



= (^»o + Hi)exp(z(5j), 



(95) 



where ~ Hi(z,i). Substituting these perturbed solu- 
tions in Eqs. (PT|) and ([5^ . and for small perturbations 
retaining only the linear terms in we get 

I 9tHi+5fHi-2AH4o(Hi+HJ) 

+ <7i2HioH2o(H2 + H*) = 0, (96) 

i 9tH2 + Aaf H2 - 2AylN4o(H2 + H*) 

+ 5i2HioH2o(H2 + H*) -0. (97) 

We consider the complex plane- wave perturbation 

Hj(z, t) = A^l cos{Kt - nz) + lAa sin(/vt - f7z) (98) 

with i = 1,2, where yl^i and Ai2 are the amplitudes for 
the real and imaginary parts, respectively, and K and f2 
are frequency and wave numbers. 

Substituting Eq. in Eqs. ([M]) and and sepa- 
rating the real and imaginary parts we get 



-AiiK = Ai2n\ 

-A12K = Aif^'- 2512^10^20^21 +4v4HtoAi 



(99) 
(100) 



for 1 = 1, and 



-A21K = A22n^X, 



(101) 



-^22^^ = ^2if^^A-25i2HioH2o-4ii+4^AH^o^2i, 



(102) 



for i = 2. Ehminating A12 between Eqs. and (jlOOp 
we get 



An[K^ - n\n^ + 4^Hto)] = -2^21512^10^201^' 



(103) 



and eliminating ^22 between (|101|) and (|102p we have 
yl2i[i^' - 172A(f}2A + 4ylAH^o)] = -2Aigi2HioH2of^'A. 



(104) 



Finally, eliminating An and ^21 from p03p and p04p 
and recalling that the density of the uniform mixture ni 
and 7i2 of the two species are given by rii = Hfg, we obtain 
the following dispersion relation 

2K = ±n[{n^ + n^X + ^An\ + ^A}?r?^ ± {(f^^ 

- ^^}? + ^An\ - ^A}?n\f + 165^2^711712}^/']^/'. 

(105) 

For stability of the plane- wave perturbation, K has to 
be real. For any f7 this happens for 

{\An\ -f ^A}?n\f > {AAnl - AAX^nlf + 165?2"i"2A, 

(106) 



or for 



4^^A71i7l2 > 5i2- 



(107) 



However, for AA^Xnin2 < 512, K can become imaginary 
and the plane-wave perturbation can grow exponentially 
with time. This is the domain of modulational instability 
of a constant-amplitude solution (uniform mixture) sig- 
nalling the possibility of coupled Fermi- Fermi bright soli- 
ton to appear [compare with inequality (fTil) of Sec. IIB 
describing stability of uniform mixture. The transforma- 
tion of the quantities in inequality (jl4p to the dimension- 
less variables of inequality (|107p can be performed with 
the definitions given after Eq. ([55)1 .] 



C. Variational Results 

Here we develop a variational localized solution to Eqs. 
((89)) and ((90)) noting that these equations can be derived 
from the Lagrangian [s^ 



L 



Ml01 + M2<^2 

X^tA/3 - 5120? 0^ 



i'l)'- (02)' -05^/3 
dz - fiiNi - fi2N2 (108) 



by demanding 5L/5(j}i ~ 6L/S(j)2 ~ 5L/Sfii = SL/5fi2 = 
0. 

To develop the variational approximation we use the 
following Gaussian ansatz [58| 



l{z) = TT" 



-1/4 



A^2Q:2 
W2 



exp(-^), (110) 
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where the variational parameters are aj, the sohtons' 
norm, and Wj width, in addition to ^j. The substitution 
of this variational ansatz in Lagrangian (jlOSp yields 



^iiVi(ai - 1) + ^2^^2(02 - 1) 



Niai N2a2 



2wl 



AalNf AXalN^ gi2NiN2aia2 
SttVSui^ Z^/inw^ yjn{'w\ + w\) 



(111) 



The first variational equations emerging from Eq. 
(fTTTj) dL/dm = dL/d^i2 = yield ai = = 1. There- 
fore the conditions ai = a2 = 1 will be substituted in the 
subsequent variational equations. The variational equa- 
tions dL/dwj = lead to 



1 



gi2N2w\ 



2NlA 



5i2iViw| 



= 0, 



= 0. 



The remaining variational equations are dL/daj 
which yield /x as a function of w-, 's, and g's: 

1 , ^/3iVl2^ , ffi2iV2 
1 



Ml 



+ 



5nwj \/Triwl + w|) ' 



^7r(u;2 + w2) 



(112) 
(113) 
0, 

(114) 
(115) 



Equations (|112p — (jllSp are the variational results which 
we shall use in our study of bright Fermi- Fermi solitons. 



D. Numerical Results 

For stationary solutions we solve time-independent 
Eqs. (|89|) and ((90|) by using an imaginary time prop- 
agation method based on the finite-difference Crank- 
Nicholson discretization scheme of time-dependent Eqs. 
(|87p and (j88p . The non-equilibrium dynamics from an 
initial stationary state is studied by solving the time- 
dependent Eqs. (j87p and ([88)1 with real time propaga- 
tion by using as initial input the solution obtained by the 
imaginary time propagation method. The reason for this 
mixed treatment is that the imaginary time propagation 
method deals with real variables only and provides very 
accurate solution of the stationary problem at low com- 
putational cost (sol . In the finite-difference discretization 
we use space step of 0.025 and time step of 0.0005. 

First we report results for stationary profiles of the lo- 
calized Fermi-Fermi solitons formed in the presence of 
attractive interspecies interaction (negative (712 )• The 
fermions form BCS state(s) which satisfy a coupled non- 
linear Schrodinger equation with repulsive (self defocus- 
ing) quintic nonlinearity. Hence fermions cannot form a 
bright soliton by itself. However, they can form a bright 
soliton in the presence of an attractive interspecies inter- 
action [3^ induced by varying an external background 
magnetic field near a Feshbach resonance [l3| ■ 
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FIG. 3: (Color online) Probability densities of the two fermion 
components from the numerical solution (labelled "num") of 
Eqs. (|89|) and (|90|) (here normalized to unity: ni{z)dz = 
1) compared with variational results (labelled "var") given by 
Eqs. (fTT2)) and (plSl) for A = 1 and (a) iVi = 10, N2 = 15 and 
<7i2 = -15, (b) Ni = 10,iV2 = 15 and 312 = -25, (c) A^i = 
im,N2 = 150 and gw = -50, and (d) Ni = 200, = 300 
and (/12 = —50. 
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FIG. 4: (Color online) Chemical potential of jth compo- 
nent obtained from the numerical solution (labelled "num") 
of Eqs. (|89|l and (|90|) and that obtained from the variational 
results (|114|l and (|115|) (labelled "var") vs. interspecies cou- 
pling pi2 for A^i = 200, N-i = 220 and A = 1. 



In Fig. [3]we present the soliton profiles of the two com- 
ponents calculated by a direct numerical solution of Eqs. 
([55)1 and (PD|) and compare them with variational results 
(|112p and (|113p . In general the numerical solutions have 
a profile distinct from a Gaussian shape of the variational 
approximation. The numerical density profile reminds of 
a square barrier. Nevertheless, the variational approx- 
imation presents a faithful average description. From 
Figs. [3] (a) and (b) we find that for a fixed A^i and N2, 
as \gi2\ is increased, the solitons become more compact 
and are better represented by variational approximation. 
From Figs. [3] (b) and (c) we see that as the number of 
fermions is increased the numerical density profiles are 
more square-barrier type than a Gaussian type. From 
Figs. [3] (c) and (d) we find that for a fixed gi2, as the 
number of atoms is reduced, the solitons become more 
compact. 

Next we illustrate how well are the variational ap- 
proximations (|114p and (|115p for the chemical poten- 
tial compared to the numerical results. In Fig. [4] 
we plot the numerically obtained chemical potential for 
Ni = 200, N2 = 220, and A = 1 for different 512 and com- 
pare with the variational result. We see that the overall 
agreement is good for all 512, although it is better for 
small I512I. 

After illustrating the soliton profiles in different states 
it is now pertinent to verify if these solitons are dynami- 
cally stable under perturbation. To this end we consider 
the typical stationary soliton of Fig. [3] (a) (obtained by 
the imaginary time propagation method) and subject it 
to the perturbation by setting (f)j{z,t) = 1.05 x (j)j{z,t) 
and observe the resultant dynamics (obtained by the real 
time propagation method). The resultant dynamics is il- 
lustrated in Fig. [5] The solitons under this perturbation 
execute some oscillation, generate some noise, neverthe- 
less propagate for as long as the numerical simulation 
was continued without being destroyed. This demon- 
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FIG. 5: (Color online) Dynamics of the probability density 
profiles of (a) the first and (b) the second Fermi solitons of 
Fig|3](a) when at i = 20 they are subject to a perturbation by 
setting (j}j{z,t) = 1.05 x (j}j{z,t). The solitons undergo stable 
propagation as long as we could continue numerical simula- 
tion. The initial soliton profile is calculated with imaginary 
time propagation algorithm and the dynamics studied with 
real time propagation algorithm. The soliton profiles are nor- 
malized to unity: nj{z,t)dz = 1. 



strates the stability of the solitons under perturbation. 
For very strong perturbation, as expected, the solitons 
are destroyed. If time covered by numerical simulation is 
too short, a unstable solution might appear to be stable. 
Thus, in numerical simulation it is important to cover 
times large compared to characteristic timescale of the 
problem, as in Fig. \5\ Also, a false stability might ap- 
pear for a small interval of time for specific space and 
time steps used in discretization. We checked the stabil- 
ity for different time steps over large intervals of time. 



V. SUMMARY 

In this paper we have obtained the phase diagram of 
a BCS superfluid Fermi-Fermi mixture of distinct mass 
fermions at zero temperature in ID, 2D, and 3D. The 
linear stability conditions relating the strength of inter- 
species Fermi-Fermi interaction with the two Fermi den- 
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sities are obtained from an energetic consideration. Two 
possible equilibrium scenarios emerge: a uniform mix- 
ture and two pure separated phases. In ID, two pure 
and separated phases appear for small fcrmion densities; 
for large densities appears the uniform mixture from an 
energetic consideration as shown in Fig. [TJ In 3D, the 
opposite happens. In addition, in 3D, a mixed and a pure 
phase can appear. In 2D, the conditions for uniform mix- 
ture and separated phases do not put any restriction on 
the fermion densities but only on the interspecies Fermi- 
Fermi interaction. 

In 3D, the uniform mixture is unstable against small 
fluctuations for large Fermi densities for a fixed 17^2 ■ For 
a positive 512 it should show partial demixing and for a 
negative gi2 it may undergo collapse. In ID, the uniform 
mixture is unstable against small fluctuations for small 
Fermi densities for a fixed 17^2 • For a positive gi2 it should 
show partial demixing and for a negative 1712 it should 
form bright solitons. Hence this mixture is of special in- 
terest for a negative 1712. This is the domain of soliton 
formation by modulational instability of the uniform mix- 
ture. To study the modulational instability and soliton 
formation in the mixture we derive a set of coupled non- 



linear equations derived as the Euler-Lagrange equation 
employing the Lagrangian density of the mixture. The 
condition of modulational instability so obtained is con- 
sistent with that of stability of uniform mixture obtained 
from an energetic consideration. In addition, we solve the 
ID dynamical equations numerically and variationally to 
study the density and chemical potential of the solitons. 
The variational result is found to be in good agreement 
with the numerical solution. We also established numer- 
ically the dynamical stability of the Fermi-Fermi solitons 
by subjecting them to a perturbation by multiplying the 
wave-function profiles by 1.05. The system is then found 
to propagate over a very long period of time without be- 
ing destroyed, which demonstrated the stability of the 
solitons. 
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